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Abstract 

We present exact results on the thermal conductivity of the one-dimensional spin- 1/2 
XXZ model in the massive antiferromagnetic and ferromagnetic regimes. The thermal 
Drude weight is calculated by a lattice path integral formulation. Numerical results for 
wide ranges of temperature and anisotropy as well as analytical results in the low and 
high temperature limits are presented. At finite temperature, the thermal Drude weight 
is finite and hence there is non-dissipative thermal transport even in the massive regime. 
At low temperature, the thermal Drude weight behaves as D{T) ^ exp(— (5/T) / \/T where 
6 is the one-spinon (respectively one-magnon) excitation energy for the antiferromagnetic 
(respectively ferromagnetic) regime. 
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1 Introduction 

Recently, transport properties of low-dimensional strongly correlated quantum systems have 
been extensively studied from both theoretical and experimental sides (see for example [1] and 
reference therein). Among them anomalously enhanced thermal conductivity [2-7] and bal- 
listic spin transport [8, 9] have been reported in experiments on one- or quasi one-dimensional 
materials with weak interchain interactions. Theoretically, the existence of such anomalous 
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properties has also been pointed out especially in quantum integrable systems. One of the 
criteria for anomalous transport is the existence of a non-zero Drude weight. 

In particular for the spin-1/2 XX Z chain, the spin transport was investigated analytically 
[10] by use of a method developed in [11]. In the massless regime, the Drude weight is non- 
zero at arbitrary temperatures resulting into non-dissipative (ballistic) spin transport. On 
the contrary, in the massive regime without magnetic field, results on the basis of [10] and 
[11] indicate a zero Drude weight for any temperature implying that the spin transport is 
dissipative. The validity, however, of these results is still controversial due to the complexity 
of the analysis at finite temperatures and debated in [12-14]. For the heat conduction, the 
thermal Drude weight in the critical regime was more recently calculated by the Bcthe ansatz 
technique in [15] which shows that the thermal transport is non-dissipative at any finite 
temperature. The study was extended by field theoretical and numerical approaches to more 
general models including non-integrable systems, and the existence of a nonzero Drude weight 
was actively discussed [16-22]. 

In this paper we discuss the thermal conductivity in the massive regime of the XXZ 
model by use of the approach developed in our previous work [15]. In the case that the 
response of a physical system to a perturbation is proportional to the force, the quantitative 
relation is obtained from linear response theory. Transport coefficients are universally given 
by the Kubo formulae [23, 24] in terms of correlation functions in thermal equilibrium without 
perturbation. In this way the thermal conductivity k{uj) relating the thermal current J7th to 
the temperature gradient, J7th = kVT, is given by the correlation function of the thermal 
current operator Jxh'- 

<^) = 7f^ die— V(*), 0(0 = / dr(Jth(-i - ir) Jth) (1-1) 
J- Jo Jo 

where /? is the reciprocal temperature; (3 = 1/T and (• • • ) denotes the thermal expectation 
value per site. Note that here we do not take into account thermomagnetic effects (cf. Sec. 5). 
The real part of eq.(l.l) reduces to 

Re k{uj) = TrDth{T)d{L0) + Kreg(w) (1.2) 

where Z)th(^^) is the thermal Drude weight^ given by 

In general, the evaluation of the above correlation function is very difficult. Fortunately, 
as already pointed out by Zotos et al. [25] , the thermal current operator of certain integrable 

systems commutes with the Hamiltonian. Indeed, the thermal current can be identified as 
one of the infinitely many conserved quantities underlying the integrability. In this fortuitous 
situation the thermal current correlations are not time dependent, they reduce to ordinary 
static correlations. Hence we find 

Ah(T) = (3^{Jl). (1.4) 

^In the previous work [15] the thermal Drude weight (written by k) is defined as k = 7r-Dth(r)- 
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Consequently, this quantity is exactly calculable within a lattice path integral formulation 
as described in [15]. At any finite temperature, the thermal Drude weight is finite implying 
non-dissipative thermal transport even in the massive regimes. At low temperature, due to 
the energy gap of the elementary excitations, the thermal Drude weight decays exponentially 
with decrease of temperature (cf. eq.(4.11) in [15] for the massless regime). 

This paper is organized as follows. In Sec. 2 we briefly review the relation between the 
thermal current operator and conserved quantities. In Sec. 3 we consider the thermal Drude 
weight at finite temperatures and discuss the results. In Sec. 4 analytical results in the low- 
and high-temperature limits are presented. Sec.5 is devoted to a summary of the present 
work and an outlook on the case of finite magnetic field. 

2 Ideal thermal current 

Let us consider the XX Z model on a periodic chain with L sites: 

^ = hkk+i hkk+i = J jf^fc (^fe+i + (^k+i^k + § {(^M+i - 1) I (2-1) 

where = ((T^±icr^)/2 and o-f. ,o'^, erf. denote the Pauli matrices acting on the kt\i space. The 
transfer integral J together with the anisotropy parameter A determine the physical nature 
of the system. In this paper we focus our attention on the massive regime A > 1, where 
correlation functions decay exponentially at zero temperature. Here J > (respectively 
J < 0) corresponds to the antifcrromagnetic (respectively ferromagnetic) model. In this 
regime it is convenient to introduce the parameter 7 instead of A: 

A = cosh 7 7GR>o. (2.2) 

Our first aim is to determine the local energy current operator j^. To achieve this we 
relate the time derivative of the local Hamiltonian to the (discrete) divergence of the thermal 
current via the continuity equation h = —divj^. Since h is written as the commutator with 
the Hamiltonian, we obtain 

hkk+1 = m, hkk+i{t)] = -{jf+i{t) - jfm. (2.3) 

Obviously the local energy current defined by 

Jk =^[hk-ik,hkk+i\ (2.4) 

satisfies the last relation in (2.3). For zero magnetic field the energy current operator Je = 
Xlfc=iife is equivalent to the thermal current operator J7e = jTth (cf. Sec.5 for non-zero 
magnetic field). Explicitly it reads 

L 

Jti. = -i^' Yl Hi4-i'^k+i - ^fe+i^fe-i) - ^(^1-1 + ^k+2)i4'^k+i - ^fc+i^fe )} • (2-5) 
fc=l 

As already shown by Zotos et.al in [25], this thermal current operator is a conserved quantity. 
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To show this from the underlying integrabihty, we consider the six vertex model which 
is the classical counterpart of the XXZ chain. There are six spin configurations carrying 
non-zero Boltzmann weights. This corresponds to six non-zero elements of the i?-matrix: 

Rlliv) = RM{v) = 1 Rll{v) = Rll{v) = Rll^y) = Rll^y) = (2.6) 

where [v] is an abbreviation for [v] = sinh(7f/2) and the meaning of the indices is exactly 
as in ref.[15]. The original quantum spin chain is connected with the classical model by a 
relation of the Hamiltonian H and the row-to-row transfer matrix T[v) = T^aYlk=i Rak{v) 



d 



dv 



2Jsinh7 



H = A—\uT{v) A = (2.7) 



i;=0 



where T{v) is a commuting family with respect to different parameters: [r{v),T[v')] = 0. 
Due to the commutativity, the transfer matrix is a generator of conserved currents J^"^"*: 

jin) ^ in-iM£))ninj^(^) D = —. (2.8) 

v=o dv 

Note that J^^^ corresponds to the Hamiltonian (2.1). For n = 2, we directly find 

L 

= i^2 ^ |4'^^^(o) - R'i,^M + [4-iik(0), i^fcit+i(0)]} R{v) := PRiv) (2.9) 
fe=i 

where P denotes the permutation operator. Thanks to the unitarity R{v)R{—v) = 1, the first 
two terms in the above equation cancel! Using the identity AR'j^j^_^^{0) = hkk+i and eq.(2.4), 
we see that J^"^^ coincides with the thermal current Jth- Due to this and the fact that the 
thermomagnetic power is always zero for zero magnetic field (being equivalent to half-filling), 
the thermal Drude weight Df^{T) is given by (1.4) (cf. (5.2) in Sec. 5). 

Here we only used the difference property and the unitarity of the i2-matrix to show the 
conservation law of the thermal current operator. These properties hold not only for the 
XXZ chain but for any integrable system whose Hamiltonian is defined as in (2.7), hence 
the thermal current operator for such an integrable system also satisfies the conservation 
law^. 



3 Thermal conductivity 

In order to obtain the thermal Drude weight (1.4), we have to evaluate the expectation value 
of the square of the thermal current operator. In fact, we are able to derive explicit rcsTilts for 
the generating function of the expectation values of any power of any conserved current. Let 
us introduce the following extended Hamiltonian H including the conserved currents (2.8) as 
a perturbation; Ti. := TC — rA,„j7*^"'' (throughout this paper we set A„ <C 1). Introducing the 
partition function 2'(A„) = Tre~'^^, we easily see that the autocorrelations of the conserved 
^One of the exceptions is the Hubbard model for which the Ji- matrix does not have the difference property. 
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quantities can be calculated by taking the second logarithmic derivative with respect to the 
variable An: 



(3.1) 



A„=0 



To evaluate Z(Xn) explicitly, we follow a procedure developed in the previous work [15]. 
Taking into account the relation (2.8), we express the partition function Z(A„) in terms of 
the row-to-row transfer matrices T{v). 



Z{Xn)= lim Trexp 

AT— »oo 



N 



nor'' i[T{uj) 



= Trexp 



N 



lim ^{lnT(^,)-lnr(0)} 



(3.2) 



where a sequence of N numbers ui, . . . , uat should be chosen such that 

N 

N^oo 



lirn^ E {In^(^i) - l^^(O)} = {-/^ + Ani"-'(^i5)"-'} lnr(t;)|^^„. (3.3) 



Applying a lattice path integral formulation, wc introduce the quantum transfer matrix 
(QTM) in the imaginary time direction [26-30]. In this formalism the partition function 
Z{Xn) and the quantity (3.1) in the thermodynamic limit L — > oo can be expressed as the 
largest eigenvalue of the QTM A: 



lim y In Z{\n) = In A ( J^")^^ - {J^'^^f 



L— »oo L 

The integral expression for A is given by^ 



d 

dXr. 



In A. 



(3.4) 



\nA = {-P + XniADr-'}£{v)l^^+ ' K{v)\u[^{v)%{v)\^ 



2iT ^ 2 cosh k'l 

fc=— oo 



(3.5) 



where £^(0) is the ground state energy of the antiferromagnetic system. Explicitly E{v) reads 

2Jsinh2 7 



E{y) = J K(v — x)e{x)dx e(i 



cos 7W — cosh 7 ' 



(3.6) 



The function 2l(v) := 1 -|- a{v) and 2l(t;) := 1 + a{v) are determined from the following set of 
non- linear integral equations (NLIEs): 

lna{v) = {-P + XniADY'-^} e{v) + k* In^iv) - k * ln^{v + 2i - ie) 

lna(t;) = {- /3 + Xn{AD)"-^] e{v) + K*ln^{v) - K*ln^{v -2i + ie). (3.7) 



For convenience, we change here the spectral parameter v — > iv. 
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Here e is an infinitesimally small number and the symbol * denotes the convolution f*g{v) = 
I-'k/'y >^(^ ~ x)g{x)dx. The kernel k{v) and the function e{v) are given by 

k(v) = -^ Y ^ e-'*^^^ e(v) = 2TrAK(v) = 2J V e''^^^^ (3.8) 

^ ^ 27r ^ 2coshA;7 ^ ^ ^ ^ ^ 2coshA;7 ^ ^ 

fc = — OO K= — CO 

From the above NLIEs (3.5) and (3.7), or just symmetry arguments, we find that expectation 
values of the conserved quantities (v7^^™^)m>i are always zero: (J^^'"^)m>i = 0. Therefore, 
due to the relations (1.4) and (3.4) (note that J^'^^ = Jth), the thermal Drude weight Ah(r) 
is given by 

2 



Ah(r) = /32(j(2)2)=/32(^^j InA. (3.9) 



We would like to remark on the structure of the NLIEs. Our NLIEs (3.7) are consistent 
with those for A„ = in [29], and may be obtained from those by the replacement of the 
driving term 

-f3e{v) + XniADr-'} e{v) (3.10) 

reflecting the structure of the general Hamiltonian Ti. and the generating function (2.8) (or 
(3.3)). Employing the same analogy, we can derive an alternative expression based on the 
thermodynamic Bethe ansatz (TBA) [31-33]. The resultant TEA equation consist of infinitely 
many NLIEs. In contrast to the TBA method, in our approach the thermal quantities are 
determined from only two NLIEs, which allows us to evaluate physical quantities numerically 
with quite high accuracy. 

In Fig.l, the temperature dependence of the thermal Drude weight Di]^{T) is depicted 
for various anisotropy parameters. We find that at low temperatures and A > 1, the weight 
I?th(2^) decays exponentially with decreasing temperature: 

Ah(r) ~ exp (-^] for T < 1. (3.11) 



where 5 is the energy gap of the one-spinon^ (respectively one-magnon) excitation in the 
antiferromagnetic (respectively ferromagnetic) regime (see the next section for details). On 
the other hand at high temperatures, £)th(T') behaves as 

Ah(T)~^ forr>l. (3.12) 

One finds that Dt\i{T) has a finite temperature maximum and the corresponding temperature 
To shifts to higher values with increasing interaction strength. In the Ising limit A — > oo, the 
temperature Tq moves to infinity for fixed J or the height of the peak goes to zero for fixed 
J A and then the thermal Drude weight Ah(r) is always zero at arbitrary temperatures. 

In Fig. 2, we show the ratio of the thermal Drude weight and the specific heat Dt\i{T)/C{T) 
as a function of the temperature T for various anisotropy parameters. Due to the mass-gap, 
this ratio goes to zero in the low-temperature limit (c.f. eq. (4.12) in [15] for the massless 
regime). In the next section the aforementioned high- and low-temperature asymptotics are 
calculated analytically. 

''Note that the excitation gap (spectral gap) actually appearing in the energy spectrum is the two-spinon 
gap 2(5. 
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4 Low- and high-temperature Hmit 



In this section we consider the low- and high-temperature behavior of the thermal Drude 
weight for both ferromagnetic (J < 0) and antiferromagnetic (J > 0) regimes. In general, 
neither the NLIEs (3.7) nor the TBA equations, can be solved analytically. However, in the 
low- and high-tcmpcraturc limits, the equations simplify especially for the gapped regimes. 
Hence we can exactly evaluate the asymptotic behavior. 



4.1 Low-temperature limit 

A. Antiferromagnetic regime ( J > 0) 

Let us consider the low-temperaturc asymptotics in the antiferromagnetic regime (J > 0). 
From the NLIEs (3.7) at low temperatures /3 3> 1, the auxiliary functions a{v) and a{v) 
reduce to 

a{v) = a{v) =exp[{-p + Xn{AD)''-'^}e{v)] for /? » 1. (4.1) 

The sub-leading terms are exponentially smaller in T than the above expression. Substituting 
(4.1) into the relation (3.5) and using (3.4), we have 

f'~ 

(j(n)2^ - = 2 r K{v) {{AD)''-^e{v)}\M-P<v)]dv. (4.2) 

7 

Applying the steepest descent method, we arrive at 

(jW2) _ (^j{n)^2 _ 2(-l)"-i(2Jsinh7) — ^ -2J«o^Binh^ 

A/-27ra2 

X \aoal + «2«Li + 2ao(«^. + «„-ia.+i) ^3^ I j 
[ AJa2 smh 7 J 

where 

«=— 00 

Thus by setting n = 2 and using (3.9), we obtain the low-temperature asymptotics of the 
Drude weight Ah(T) 

^ -2(2Jsinh7)taoa2„.g^i^ r^_| ^ ^ X | _ 

Here we have used the fact a2m-i = when m > 1. Note that the exponent — 2Jaosinh7 
in (4.5) is nothing but the one-spinon excitation energy of the model. This result indicates 

that the spinon excitation mainly contributes to the low-temperature heat conduction in 
the antiferromagnetic regime. Comparing the low-temperature limit for the specific heat 
C7(T) = /32((j(i)2) _ (^(i))2) 

^ 2(2Jsinh7)tag„.g:£^ f^.a ^ 3 + O hh)] (4.6) 

V-27ra2 I 4jQ;osinh7 ^ 
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we have 

= \ -T + 0{T^). 4.7 

In the inset of Fig. 2(a), we show these results for various anisotropy parameters. 
B. Ferromagnetic regime (J < 0) 

Next wc analyze the low-temperature asymptotics in the ferromagnetic regime (J < 0). 
In this regime the auxiliary functions a{v) and a{v) are no longer small at low-temperature, 
hence the analysis of the NLIEs becomes difficult. To deal with this situation, we introduce 
the alternative equations by taking the reciprocals [29] : 

b{v) = biv) 



a{v) a{v) 
^{v) = 1 + b{v) ^{v) = 1 + b{v). (4.8) 

Applying the Fourier transform, we derive the NLIEs in terms of these functions: 

lnb{v) = {P- XniADY'-^} e{v) + k * ln*B(t;) - k * ln^(t> + 2i - ie) 

lnb{v) = {/3 -A„(^L»)"-^} £*(?;) + K* In ^(v) -K* In <B(t;-2i + ie) (4.9) 

where the driving terms £{v) and e*{v) are given by 

~. s cosh ^^(1 — i?;) , . cosh 77(1+ if) 

^W = -"^ 2si.hi(l-„) ^W-^^ 2s.nhi(l + J) 

and the kernel k{v) is defined as 

^(^) = _^y£:!l^^!MiM. (4.11) 

^ ' 27r^ sinhA;7 ^ ^ 

k=i ' 

Accordingly, the largest eigenvalue A is written as 

TT TV 

lnA= K{v)ln^{v)dv+ [ K*{v)lnm{v)dv 

Applying an iteration procedure to (4.9), we obtain the low-temperature behavior of the 
auxiliary functions: 

TT TT 

In b{v) = g{v) + I'k{v- x)e^(^)dx - / K{v + 2i- ie)e^*('')da; 

1 T 
/. — /» — 

\u.b{y) = g* {v) + ' li{v-x)e^'^^^dx- j k(v - 21 + ie)e^^^^ dx (4.13) 
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where 

g{v) = {13- Xn{ADr-^}e{v) g*{v) = {(3- Xn{ADr-^} T {v) . (4.14) 

Identifying the poles of 'k{v) a,t v = ib2i and applying Cauchy's theorem, the dominant 
contributions of the integrals for T <C 1 are evaluated 

/— /•— — -7(i>— a;+ie) 
K(?;-x + 2i-ie)ln(l + e^*(^))dx- --^ r ^4- -e»*(^)dx = -e^*(^\ 
-I 27r J_E. 2 sinh ^'yyu — x + ie) 

(4.15) 

Using this and neglecting the subdominant terms, we obtain 

ln«8(u) ~ b(t;) - e^'(^)(l + e^'*^^)) = e^'(^) + e^'^^^+s*^''). (4.16) 
For ln©(w) a similar equation is valid. Substituting the results into (4.12), we have 

J -TL J-K J-2L 

We can calculate the first (respectively second) term in (4.17) by shifting the contour to — ioo 
(respectively ioo). The third integral is evaluated by a saddle point integration. 

Using the relation (3.4), we obtain the low-temperature asymptotics of the conserved 
quantities: 

1 4n — 1 

V— 27rp2 



where 



—4 J/32 sinh 7 



Setting n = 2, we derive the low-temperature asymptotics of the thermal Drude weight 
Ah(T): 



2./(A-l) 



o.,(r)=(z2:Qli^^iI^^{r-i + o(Ti)}. (4.20) 



2n 

The exponent 2|J|(A — 1) in (4.20) agrees with the one-magnon gap. This result coincides 
with that from numerical diagonalization of the Hamiltonian for size L = 18 [20]. For the 
specific heat C{T) = /^^((^(i)^) - (j(i))2) we obtain 

,o Jsinh7 , '-i , ^„ 2J(A-1) ^ 

J sinh 7^6^^ (-2J)2(A- l)2e^^ ( 3 3 _ 1 

(4.21) 

Here the exponent \ J\ sinh7 in the specific heat (4.21) is the one-spinon excitation gap in the 
ferromagnetic regime. We observe a crossover behavior from dominant one-magnon excitation 



9 



(A > 5/3) and dominant one-spinon excitation (A < 5/3) [33,34]. In contrast to this, a 
crossover does not take place in the low-temperature thermal Drude weight, which implies 
that only the magnon excitation contributes to the low-temperature heat conduction. This 
result yields a different behavior of the ratio of the thermal Drude weight and the specific 
heat (4.22) for small and large values of A, respectively 



Ah(T) 
C{T) 



-2JT + 0{T'^) 
(-2J)i(A- l)e' 



for A < -, 



27rJ2(A + l) 



|Ti+0(|Ti)}, forA> 



(4.22) 



The ratio of -Dth(^) and C(T) depends linearly on the temperature and is independent of the 
anisotropy for A < 5/3. On the other hand, for A > 5/3 the ratio decays exponentially with 
temperature and explicitly depends on A. This behavior is clearly observed in Fig. 2(b). 



4.2 High temperature limit 

Here we analyze the high-temperature behavior. In this limit the NLIEs (3.7) linearize and 
hence can be analytically solved. By use of identities like 



— lno= -— ln2t 
d/3 a dp 



52 



21 92 



In 21 



(4.23) 

(4.24) 



and the limiting behavior a = 1, 21 = 2 for /? = we find the linearized integral equations 



^ In^iv) = ~e{v) + ln2l(t;) " ^'^ * ^ ln2l(i; + 2i - ie) 

ln2l(^;) = -^e{v) + * ^ ln2l(i;) - ^« * ^ ln2l(i; - 2i + ie) 



d_ 

dp 



(4.25) 



and integral equations linear in [d/dXf' In 21 (with A := An) 



d_ 



— In2l(i;)j — j ln^(v)--K*( — \ ln2l(i; + 2i - ie) 



A)^n^(.)=(|-ln2l(/' 



(4.26) 



The current correlator (3.1) is found from (3.4) 



dv 



A=0 



dv + c.c. 



(4.27) 



A=0 
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where in the last hne we have used the dressed function formahsm. With the high temperature 
asymptotics a{v),a{v) = 1 for /3 = we find 



dX 



In A 



A=0 



8ttA J_n df3 \ dX J 



dv + c.c. 



(4.28) 



A=0 



The integrands in the above equation are found analytically. First, from (4.25) we obtain 



da{v) 



dp 



A=0 



1 



J sinh^ 7 

2 sin ^{v + i) \ sin ^{v + 3i) sin ^{v — i] 



(4.29) 



and a is the complex conjugate. 

Similar to the reasoning at the beginning of this section we obtain linear integral equations 
for the derivatives with respect to A 



d (AD)"'~^ 1 d 

- In2l(.) = '-^e{v) + * - ln2t(.J - - 



19- 
K* — ln2l(w + 2i) 



dX 



d — (AD)'^~^ Id— Id 

_ in2t(^;) = ^ e{v) + * ^ ^^^i^) - * g^^^^i^ - '^^) ■ 



Hence we find in the high temperature limit 



d_ 
dX 



a{v) 



A=0 



A=0 



-{ADr-^±a{v) 



A=0 



(4.30) 

(4.31) 
(4.32) 



A=0 



By use of these explicit expressions for n = 2 we obtain the high temperature limit of Di^{T) 



Ah(T) 



J^(2 + cosh27) 1 

2 r2 



+ o 



1 

T3 



(4.33) 



which is consistent with the result in [15]^ after changing the parameter 7 — > 17. 



5 Summary and discussion 

In this paper we have discussed the thermal transport properties in the massive regime for 
the spin-1/2 XX Z chain. The thermal current operator is expressed as a conserved quantity 
resulting into an anomalous thermal transport. The thermal Drude weight at finite tem- 
peratures was calculated by a path integral formulation. Due to finite temperature effects, 
non-dissipative thermal transport was observed even in the massive regime of the model. At 
low-temperatures (T <C 1), Dih{T) can be written in a universal form: Dth(?') ~ e~^^'^ /VT, 
where S is the one-spinon (respectively one-magnon) gap for the antiferromagnetic (respec- 
tively ferromagnetic) regime. 

Finally, we would like to remark some generalizations of our results, (i) As mentioned 
in Sec. 2, the conservation law of the thermal current is not limited to the present model. 

^There is a misprint in (4.18) in [15]: on the r.h.s of (4.18) a factor tt is missing. 
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Therefore the present approach is directly apphcable to more general models such as the 
XYZ and the integrable higher-spin or higher-rank chains, (ii) It is quite interesting to 
consider the system in an external magnetic field h: H = TC + h (t| /2. From the definition 
(2.4), it is observed that the energy current J^e includes the spin current J'g proportional to 
h. Namely 

L 

Je = Jth + hj, J, = iJ ^(<t^7+i - 4+1^7) (^-1) 

j=i 

where the explicit form of Jth is given by (2.5). Despite the fact that the thermal current 

is still conserved \H, i7th] = 0, the energy current is no longer a conserved quantity because 
[H^ i7s] 7^ 0. To evaluate the thermal conductivity at /i > 0, we should take carefully into 
account thermomagnetic effects, which does not matter when /i = (being equivalent to 
half-filling) because {JsJtii) = 0. The correct thermal Drude weight at ^ > is determined 
from linear-response theory (see [24] for example): 

Ah(r, h) = (3^{jl)n - ^S{T, hfD,{T, h) S{T, h) = P'^-^^ (5-2) 

where {■ ■ ■)h denotes the thermal expectation value per site for the system H; S{T, h) is the 
thermomagnetic power and Ds{T,h) is the Drude weight of the spin stiffness a{u): 



Re £7(0;) = TtDs{T, h)5{uj) + C7reg(w) 



D,{T,h) = ^{ 



{-K)h-2 V 



P / 171 m n 



\{n\Js\'m)\ 



2 



Pn = — (5.3) 



Here, K and E^^ denote the kinetic energy and the energy eigenvalues of the Hamiltonian H, 
respectively. The first term in (5.2) can be evaluated by an extension of the present method. 
The Drude weight Ds(T,h) in the second term may be calculated by the approach given in 
[11]. In fact in the massless regime (A < 0) with zero external field, Ds(T, 0) has been already 
calculated by Zotos [10]. As mentioned in the preceding section, however, the validity of the 
resultant Drude weight is presently debated (see [1, 12-14] for example). In this respect, a 
rigorous study of the Drude weight together with the thermomagnetic power is highly desired 
to investigate the thermal transport in finite magnetic fields. 



Acknowledgments 

We would like to thank M. Shiroishi, M. Takahashi and J. Takeya for stimulating discussions. 
The authors acknowledge financial support by the Deutsche Forschungsgemeinschaft under 
grant No. Kl 645/4-1 and SP1073. KS is supported in part by the JSPS research fellowships 
for young scientists. 



12 



References 

[I] X. Zotos and P. Prelovsek, cond-mat/0304630. 

[2] K. Kudo, S. Ishikawa, T. Noji, T. Adachi, Y. Koike, K. Maki, S. Tsuji and K. Kumagai, 
J. Low. Temp. Phys. 117, 1689 (1999). 

[3] K. Kudo, S. Ishikawa, T. Noji, T. Adachi, Y. Koike, K. Maki, S. Tsuji and K. Kumagai, 
J. Phys. Soc. Jpn. 70 437 (2001). 

[4] A. V. Sologubenko, K. Gianno, H. R. Ott, U. Ammerahl and A. Revcolevschi, Phys. Rev. 
Lett. 84, 2714 (2000). 

[5] A. V. Sologubenko, K. Gianno, H. R. Ott, A. Vietkine and A. Revcolevschi, Phys. Rev. 
B 64, 054412 (2001). 

[6] C. Hess, C. Baumann, U. Ammerahl, B. Biichner, F. Heidrich-Meisner, W. Brenig and 
A. Revcolevschi, Phys. Rev. B 64, 184305 (2001). 

[7] A. V. Sologubenko, K. Gianno, H. R. Ott, U. Ammerahl and A. Revcolevschi, Phys. Rev. 
Lett. 84, 2714 (2000). 

[8] M. Takigawa, N. Motoyama, H. Eisaki and S. Uchida, Phys. Rev. Lett. 76, 4612 (1996). 
[9] K.R. Thurber, A.W. Hunt, T. Imai and F.C. Chou, Phys. Rev. Lett. 87, 247202 (2001). 
[10] X. Zotos, Phys. Rev. Lett. 82, 1764 (1999). 

[II] S. Fujimoto and N. Kawakami, J. Phys. A 31, 465 (1998). 

[12] B. N. Narozhny, A. J. Millis and N. Andrei, Phys. Rev. B 58, R2921 (1998). 

[13] J. V. Alvarez and C. Gros, Phys. Rev. Lett. 88, 077203 (2002). 

[14] S. Fujimoto and N. Kawakami, Phys. Rev. Lett. 90, 197202 (2003). 

[15] A. Kliimper and K. Sakai, J. Phys. A: Math. Gen. 35 2173 (2002). 

[16] J. V. Alvarez and C. Gros, Phys. Rev. Lett. 89, 156603 (2002). 

[17] K. Saito and S. Miyashita, J. Phys. Soc. Jpn. 71, 2485 (2002). 

[18] K. Saito, Phys. Rev. B 67, 064410 (2003). 

[19] K. Saito, Phys. Rev. B 67, 164410 (2003). 

[20] F. Heidrich-Meisner, A. Honecker, D.C. Cabra and W. Brenig, Phys. Rev. B 66, 
140406(R) (2002). 

[21] F. Heidrich-Meisner, A. Honecker, D.C. Cabra and W. Brenig, cond-mat/0304595. 
[22] E. Orignac, R. Chitra and R.Citro, Phys. Rev. B 67, 134426 (2003). 



13 



[23] R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957). 

[24] G. D. Mahan, Many-Particle Physics (New York: Plenum Press). 

[25] X. Zotos, F. Naef and P. Prelovsek, Phys. Rev. B 55, 11029 (1997). 

[26] M. Suzuki, Phys. Rev. B 31 (1985) 2957. 

[27] M. Suzuki and M. Inoue, Prog. Theor. Phys. 78 (1987) 787. 

[28] A. Kliimper, Ann. Physik 1, 540 (1992). 

[29] A. Kliimper, Z. Phys. B 91, 507 (1993). 

[30] C. Destri and H.J. de Vega, Phys. Rev. Lett. 69, 2313 (1992). 

[31] M. Gaudin, Phys. Rev. Lett. 26, 1301 (1971). 

[32] M. Takahashi, Prog. Theor. Phys. 46, 401 (1971). 

[33] M. Takahashi, Thermodynamics of One-Dimensional Solvable Models, (Gambridge: 
Cambridge University Press). 

[34] M. Takahashi, Prog. Theor. Phys. 50, 1519 (1973). 



14 



1.4 



1.2 



0.8 



0.6 



0.4 



0.2 



1 1 

- (a) J>0 


1 

^=1.0 




A^\.5 








A=3.0 




\ A=A.Q ---- 


- / 




r 1 1 



0.5 



T/J 



0.35 



0.3 



0.25 



0.2 - 



q 0.15 



0.1 - 



0.05 



0.5 



T/J 



1.5 



(b) /<0 




A^\.Q5 




A-\.\Q 






^-1.20 






^ A^IAO 






A-im 






^-4.00 ---- - 











1.5 



Figure 1: Illustration of numerical results for the thermal Drude weight Di^{T) (in units of 
(JA)^) in the antiferromagnetic J > (a) and ferromagnetic J < (b) regimes as a function 
of temperature T (in units of J) for various anisotropies. 
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Figure 2: Illustration of the ratio of the thermal Drude weight and the specific heat 
Dth{T)/C{T) (in units of J^) in the antiferromagnetic J > (a) and ferromagnetic J < 
(b) regimes as a function of temperature T (in units of J) for various anisotropics. In the 
inset numerical data at low temperatures are shown. The linear lines in the inset of Fig.2(a) 
are the analytical results for the low-temperature asymptotics. 



16 



